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Abstract. The gravitino dark matter with neutralino NLSP hypothesis is investigated in 
the framework of NMSSM. We have considered both the thermal and non-thermal gravitino 
production mechanisms, and we have taken into account all the collider and cosmological 
constraints. The maximum allowed reheating temperature after inflation, as well as the 
maximum allowed gravitino mass are determined. 



1. Introduction 

There is accumulated evidence both from astrophysics and cosmology that about 1/4 of the 
energy budget of the universe consists of so called dark matter pQ, namely a component which 
is non-relativistic and neither feels the electromagnetic nor the strong interaction. It is known 
that in the standard model of particle physics there is not a suitable dark matter candidate. 
Supersymmetry, a very well-motivated idea, is perhaps the most popular way beyond the 
standard model, and it provides us with an ideal candidate for playing the role of cold dark 
matter in the universe. 

Out of the several dark matter candidates, certainly the most well-studied case is the one of 
the neutralino. However, the gravitino is also a very compelling candidate for several reasons. 
Its interactions are completely fixed by the supergravity Lagrangian, its mass is directly 
related to the scale at which supersymmetry is broken, and finally the reheating temperature 
after inflation Tr, which is essential for the precise baryon asymmetry generation mechanism, 
plays a certain role in gravitino cosmology. 

The MSSM [2] is the most economical supersymmetric extension of the standarm model. 
However, it suffers from the so-called fi problem [3J. The natural values for the higgs/higgsino 
mass parameter [i is either or the Planck mass, while phenomenology requires the fi should be 
at the electroweak scale. The simplest supersymmetric model that solves the n problem is the 
NMSSM [4], where one more singlet chiral superfield is introduced. 

It is therefore interesting to try to see whether a particular cosmological scenario, with 
gravitino LSP and neutralino NLSP in the framework of the NMSSM, is a viable one. 



2. The theoretical framework 

2.1. Basics of NMSSM 

The particle physics model is defined by the superpotential 

W = e tj (y u H{ Qiu + Ya H\ Q j d + Y e H\ V e) - e^A S H\h{ + ^kS 3 (1) 



£ so ft = rnA Q* Q + m| u* tt + m| d* d + m| L* L + m| e* e 



as well as the soft breaking masses and couplings 

' y: ^ \- u*u + m 2 ~ 

+m 2 Hi H* Hx + mf^ // 2 *^ 2 + m 2 s S*S 
+eij U u Y u H{ tyu + Ad Y d H\ Qi d + A e Y e H\ V e + H.c. 

+ (-eyA A X SH\H 3 2 + i K ,4 K 5 3 + H.c.^ 
- i (M 3 A 3 A 3 + M 2 A 2 A 2 + Mi Ai Ai + H.c. ) (2) 

When the singlet acquires a vaccum expectation value, S, we obtain an effective fi parameter, 
H e ff = XS. Imposing universality at the GUT scale, a small controllable number of free 
parameters remains, namely 

tan(3 = v u /v d , m , A ,m 1 / 2 , A, A k 
and the sign of the effective /x parameter can be chosen at will. 

Because of the extra singlet superfield, in the NMSSM there is a larger higgs sector and a 
larger neutralino sector. The neutralino mass matrix is characterized by the appearence of a fifth 
neutralino state, meaning that the composition of the lightest neutralino has an extra singlino 
contribution 

X? = iVii^ + iVi2^3° + N 13 H° + N U H° 2 + N 15 S (3) 

In the following, neutralinos with N 2 -^ > 0.9, or N 2 5 > 0.9, will be referred to as bino- or 
singlino-like, respectively. 

Furthermore, in the Higgs sector we have now two CP-odd neutral, and three CP-even neutral 
Higgses. We make the assumption that there is no CP-violation in the Higgs sector, and therefore 
the CP-even and CP-odd states do not mix. We are not interested in the CP-odd states, while 
the CP-even Higgs interaction and physical eigenstates are related by the transformation 

h° a = S ab H° b (4) 

where S is the unitary matrix that diagonalises the CP-even symmetric mass matrix, a,b = 
1,2,3, and the physical eigenstates are ordered as m h o < m h o < m h o. 

2. 2. Production of gravitinos 

After inflation gravitinos can be produced in two ways. One way to produce gravitinos is with 
scatterings from the thermal bath, and another is from the out-of-equillibrium decays of the 
NLSP, which decouple from the thermal bath before primordial Big-Bang Nucleosynthesis and 
decay after the BBN time. Thus, imposing the WMAP bounds [5] we can write for the gravitino 
abundance 

0.1097 < fl 3/2 h 2 = ^f 2 h 2 + Y^ /2 LSP h 2 < 0.1165 (5) 

where 

n NLS Ph2 = ( 
' rn NL sp 



with rriNLSP the mass of the NLSP, and 0<NLSph 2 the abundance the NLSP would have, had it 
not decayed into the gravitino. The thermal contribution is given by (approximately for a light 
gravitino, tHq <C rrig) [6] 

2.3. Neutralino life-time 

The NSLP is unstable with a lifetime that is typically larger than BBN time tsBN ~ 1 sec. 
Energetic particles produced by the NLSP decay may dissociate the background nuclei and 
significantly affect the primordial abundances of light elements. If such processes occur with 
sizable rates, the predictions of the standard BBN scenario would be altered and the success 
of the primordial nucleosynthesis would be spoiled. BBN constraints on cosmological scenarios 
with exotic long-lived particles predicted by physics beyond the Standard Model have been 
studied, and here we have used the figure 2 of Ref. [7]. Thus, we need to compute the lifetime of 
the neutralino. There are three main decay channels, namely neutralino decaying into gravitino 
and a standard model particle that can be one of the followings: Photon, Z boson and Higgs. 
The supergravity Lagrangian is known, and the Feynman rules for the gravitino interactions 
have been derived [8]. Therefore, the neutralino lifetime is given by 

r = ~ (8) 

r = r(x -> iG) + r( x -> zg) + r( x hG) (9) 

and the exact expressions for the diferent contributions can be found in [S]. 
3. Constraints and results 

We have used NMSSMTools [TO], a computer software that computes the masses of the Higgses 
and the superpartners, the couplings, and the relic density of the neutralino, for a given set of 
the free parameters. We have performed a random scan over the whole parameter space (with 
fixed /x > motivated by the muon anomalous magnetic moment), and we have selected only 
those points that satisfy i) theoretical requirements, such as neutralino LSP, correct electroweak 
symmetry breaking, absence of tachyonic masses etc., and ii) LEP bounds on the Higgs mass, 
collider bounds on SUSY particle masses, and experimental data from B-physics. Finally, for 
any given point in the cNMSSM parameter space, the neutralino lifetime is a function of the 
gravitino mass only. Imposing the BBN constraints we find the maximum allowed gravitino 
mass, and from the cold dark matter bound we can determine the maximum allowed reheating 
temperature. For all the acceptable points the lightest neutralino is either a bino or a singlino, 
and our main results are summarized in the figures below. 

For the bino case, figure 1 shows the maximum allowed reheating temperature after inflation 
versus the maximum allowed gravitino mass, both in GeV. Although it cannot be seen directly 
from the figures, the maximum possible gravitino mass in the bino case is mQ — 1 GeV, and the 
corresponding reheating temperature is Tr ~ 10 GeV. Therefore, we see that a) the gravitino 
in this scenario must be much lighter than the rest of superpatners, but nevertheless it is still in 
the gravity-mediated SUSY-breaking scheme, and b) the reheating temperature after inflation 
is not large enough for thermal leptogenesis. 

For the singlino case, we show in figure 2 the maximum allowed reheating temperature versus 
gravitino mass, both in GeV. This time the neutralino relic density is even larger than before, and 
gravitino now must be extremely light. This is due to the smallness of the coefficients N\i, N\2 
in the decay rate to gravitino and photon. For the same lifetime as before, the gravitino mass 
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Figure 1. Reheating temperature versus gravitino mass (both in GeV) for the bino case. 
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Figure 2. Reheating temperature versus gravitino mass (both in GeV) for the singlino case. 

must be several orders of magnitude lower than in the bino case. The last figure shows that 
in the singlino case the reheating temperature cannot be larger than about 200 GeV. However, 
this value is much lower than the minimum value required for the computation of the gravitino 
thermal production Tr > 1 TeV, and therefore we conclude that this scenario must be excluded. 

4. Conclusion 

In the framework of the cNMSSM, we have considered a possible cosmological scenario with the 
gravitino LSP and the neutralino NLSP. The gravitino is stable and plays the role of cold dark 
matter in the universe, while the neutralino is unstable and decays to gravitino. We have taken 
into account the relevant gravitino production mechanisms, which are i) the NLSP decay, and ii) 



scattering processes from the thermal bath. Our results can be seen in the figures. We have found 
that i) the gravitino is necessarily very light, and ii) the reheating temperature after inflation 
is two orders of magnitude lower than the temperature required for thermal leptogenesis. The 
singlino scenario must be excluded, while in the bino case it is hardly possible to have a gravitino 
in the gravity- mediated SUSY-breaking scheme. 
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